Proximity effect in ultrathin Pb/Ag multilayers within the Cooper limit 
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We report on transport and tunneling measurements performed on ultra-thin Pb/Ag (strong cou- 
pled superconductor/normal metal) multilayers evaporated by quench condensation. The critical 
temperature and energy gap of the heterostructures oscillate with addition of each layer, demonstrat- 
ing the validity of the Cooper limit model in the case of multilayers. We observe excellent agreement 
with a simple theory for samples with layer thickness larger than 30 A. Samples with single layers 
thinner than 30A deviate from the Cooper limit theory. We suggest that this is due to the "in- 
verse proximity effect" where the normal metal electrons improve screening in the superconducting 
ultrathin layer and thus enhance the critical temperature. 
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The superconducting proximity effect is a well known 
phenomenon and has drawn a lot of interest both from 
the fundamental and the practical points of view. In 
a high transmission contact between a superconductor 
and a normal metal the superconducting wave-function 
varies smoothly across the interface causing a suppres- 
sion of the pair amplitude in the superconductor and an 
enhancement of pairing on the normal side. The charac- 
teristic distance in which superconductivity "leaks" into 
the normal region is the normal state coherence length 
— \/hD/2n/cT, where D is the diffusivity, and the 
length-scale at which superconductivity is suppressed in 
the superconducting side is the superconducting coher- 
ence length, ^s- 

In the Cooper limit JJ (in which both the supercon- 
ductor and the normal metal layers are smaller than the 
characteristic coherence lengths and respectively) 
the behavior of an S-N bilayer is well described within the 
framework of the de Gennes model 0, Q • The electrons 
experience an average paring interaction (average be- 
tween the two materials) and Tc can be described within 
the BCS weak coupling form adapted to the Cooper limit 
proximity effect in the following way : 
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here cod is the debye frequency, N(0)V is the pairing in- 
teraction (N(0) being the density of states at the Fermi 
level and V the pairing interaction), and ds,N is the thick- 
ness of the superconducting and normal metal films re- 
spectively. This model assumes that the debye frequen- 
cies in the two metals are similar and it is based on the 
understanding that the electrons spend 

their time in the normal metal and 



{Nr^dN+Nsds) 
Nsds ;„ 4-\,„ 



superconductor {Nn and Ns are the density of states 
in the normal metal and in the superconductor respec- 
tively) . 



The above picture is modified in the case of strong cou- 
pled superconductors where Tc is given by the McMillan 
expression |^ : 
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here A is the dimensionless electron-phonon coupling 
coefficient and fi* is the effective coulomb repulsion. In a 
normal- metal/strong coupled superconductor within the 
Cooper limit one can expect A in equation |21 to be re- 
placed by an average electron-phonon coupling : 
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The proximity effect in the Cooper limit has been ob- 
served experimentally in numerous systems of supercon- 
ductor/normal metal bilayers. In this paper we study 
systems of multilayers of superconductor/normal metal 
(Pb/Ag) in which the total thicknesses of the metals are 
comparable to the relevant coherence lengths. The trans- 
port and tunneling experimental results confirm that the 
basic idea of the Cooper limit notion, i.e. that the elec- 
trons experience an average pairing, is valid also in a 
multilayer sample. We observe good quantitative agree- 
ment with theory as long as the thickness of the layers 
is larger than ^ 30A while thinner layers deviate from 
the theoretical predictions. We suggest that the simple 
picture is modified in ultra-thin layers due to effects of 
screening of e-e interactions (the inverse proximity effect 
as previously reported 0). 

The Pb and Ag layers were prepared by quench con- 
densation, i.e. evaporation on a cryogenically cooled 
substrate within the measurement apparatus |3- This 
allows in-situ sequential depositions under UHV condi- 
tions and simultaneous transport and tunneling measure- 
ments on a single sample. This method has several es- 
sential advantages for the study of the superconducting 
proximity effect Since Pb and Ag are immiscible, 
the alloying of these two materials is very improbable 
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FIG. 1: Resistance per square versus temperature of a Pb/Ag 
bilayer (sample #2). The thickness of Ag varies from A 
(Pb) to 14.4 A (Ag9). Each evaporation corresponds to ap- 
proximately 1.6 A. The Tc decreases monotonically as normal 
state resistance is decreased. Insert: Illustration of the sample 
geometry. 

even at room temperatures. Alloying becomes practica- 
bly negligible when the samples are quench condensed 
and are kept below T=10K throughout the entire exper- 
iment (sample growth and measurement). In addition, 
the extremely clean environment of the experiment leads 
to very high quality interfaces between the superconduct- 
ing layers and the normal metal layers. Such clean inter- 
faces provide ideal conditions for proximity effects. 

The geometry of the samples used in this study is 
shown in the insert of figure 1. We begin by evaporat- 
ing two strips of Al onto a room temperature polished 
glass substrate. These strips are used both as voltage 
leads and as the base electrode for tunneling measure- 
ments. Then we allow a native oxide barrier to grow on 
the Al for 20 minutes. We connect metallic leads to the 
substrate and place it in an evaporation chamber which 
is then pumped out and immersed in liquid helium. A 
lOA Ge adhesion layer is then quench condensed across 
the Al strips, permitting the evaporation of a continu- 
ous layer of Pb as thin as SA 0. The thicknesses are 
monitored by a calibrated quartz crystal situated in the 
chamber. Multilayers, made of base units of Pb/Ag bi- 
layers having different thicknesses, are quench condensed 
across the Al strips. This configuration allows the si- 
multaneous measurements of transport properties of the 
multilayers (thus determining the critical temperature), 
and tunneling into the multilayers (for determination of 
the energy gap). We performed these measurements at 
incremental stages of the evaporation of each of the Pb or 
Ag layers. All measurement were performed in a screened 
room using standard AC methods. 

In this paper we present data on three Pb/Ag systems 
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FIG. 2: Superconducting critical temperatures extracted from 
the R(T) measurements in the case of three different thickness 
of the Pb/Ag layer. The thicknesses of the Pb and Ag layers 
in the bilayer are 32.5 A (top graph), 14.4 A (middle graph) 
and 9.3 A (bottom graph). The solid lines are the theoretical 
fits using equations 2 and 3. 



having different base unit thickness: Sample #1 where 
dpb = dAg ~32A, sample #2 where dpb = dAg ~15A, 
and sample #3 where dpb = dAg ~9A. Each of these 
unit bilaycrs is within the cooper limit; of the Ag at 
low temperatures is of the order of 1 /im and ^5 of the 
amorphous Pb (which is in the dirty limit) at low tem- 
perature is given by ^' = VSfll where ^0 is the clean limit 
coherence length (for Pb this is 8OOA) and I is the mean 
free path which is approximately an interatomic distance, 
about 2.7A in our samples 0. These values yield a zero 
temperature coherence length, ^' « 46A, which is larger 
than each one of our Pb layers. Furthermore, as the tem- 
perature increases the coherence length increases until it 
diverges at T — Tc- 

We begin by presenting the results for the initial bi- 
layer. Figure 1 shows the resistance versus temperature 
of sample ^2 as a layer of 14.5 A thick Ag is grown in 
steps of « 1.55 A on top of a 14. 4A thick initial Pb layer. 
Such measurements are used to extract Tc by taking the 
midpoint of the transition {'Rn/2) j^. It is seen that the 
initial Pb layer has a Tc which is smaller than the bulk 
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value of 7.2K. This is typical for very thin superconduct- 
ing layers which are known to exhibit a superconductor 
- insulator transition as a function of film thickness (and 
Rs?) milllll- It has been argued that as R^, 

(the disorder) increases, the electronic screening is re- 
duced. Coulomb interactions are then enhanced leading 
to a weakening of the superconducting coupling. This 
explains the decrease of Tc in the Pb films when the 
thickness is reduced from 30A to QA as is seen in figure 
2. The enhancement of Coulomb interactions for high 
Ksq is also evident in figure 3 which depicts the dl/dV-V 
measurements (proportional to the density of states) of 
a 9A thick Pb film with increasing layers of Ag, taken at 
temperatures above Tc- It is seen that in the ultrathin 
Pb layer the density of states at the Fermi level is rela- 
tively low. The tunneling conductance exhibits a strong 
voltage (energy) dependence demonstrating that this is 
a strongly coulomb correlated system. The density of 
states dramatically increases when layers of Ag are evap- 
orated on top of the Pb, presumably due to enhancement 
of screeningand, hence, suppression of the Coulomb in- 
teractions 

As subsequent layers of Ag are added (and R^g is re- 
duced as 1/ dAg ) Tc is decreased even further due to the 
proximity effect as evident in figure 1. Since the bilayer is 
in the cooper limit, Tc is indeed expected to decrease as 
a function of the normal metal thickness. We note, how- 
ever, that this is not always the case. We have shown jij 
that when the initial Pb layer is extremely thin (thinner 
than that of sample 3) the first added Ag sub-layers cause 
Tc to increase in what seems to be an inverse proximity 
effect . This behavior was interpreted as being due to the 
fact that the Ag layers enhance screening in the system 
thus suppressing the e-e interactions which are responsi- 
ble for the Tc suppression. This trend opposes the usual 
proximity effect and becomes more important in thinner 
superconducting films. 

The evolution of Tc as more Pb and Ag layers are 
added is shown in figure 2 for the three samples. It is 
seen that Tc oscillates when we change between Pb and 
Ag depositions reaching a local minimum or maximum 
at the completion of each layer. The oscillation ampli- 
tude decreases with increasing number of bilayers and 
Tc appears to approach a value of « 2iK in all three 
samples. The oscillation behavior is also observed in the 
energy gap measurements which are illustrated in figure 
4 for sample #1. The dl/dV curves were measured at 
T==1.65K well below the Tc of the multilayer. Note that 
despite that fact that tunneling is performed into the bot- 
tom Pb layer, the results are still affected by adding the 
6th layer on top. This clearly demonstrates that within 
the cooper limit the tunneling electrons probe the entire 
sample and the superconducting parameters are deter- 
mined by the mean properties of all the layers in the 
multilayer. Figure 4 also shows the ratio as a func- 
tion of the number of layers. The superconducting gap 
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FIG. 3: dl/dV versus V (proportional to the density of states 
as a function of energy) at temperatures above Tc (2.Q4K) of 
an ultrathin Pb layer (10 A) with an increasing Ag (up to 10 
A) layer on top. 



was evaluated from a classical temperature dependent 
BCS fit. It is seen that this ratio crosses over from a 
value of 4.7 (which is close to the strong coupling limit 
of ~4.8) to 3.6 (which is close to the week coupling limit 
of ~3.5). The fact that this ratio is relatively constant 
as we add the last few layers gives us further confidence 
that all our multilayers are within the Cooper limit. 

We have attempted to fit the curves in figure 2 to the 
proximity effect modified McMillan expression of equa- 
tion [3 using an average A for the Cooper limit (equa- 
tion in the spirit of the de-Gennes model. These are 
shown in figure 2 were we used the known values for 
Pb: As = 1.55, W£) = 105X, ^* = 0.11 and taking a low 
electron-phonon coupling for Ag, \n — 0.2, without any 
fitting parameters. It is seen that the experimental re- 
sults are in excellent agreement with this model for sam- 
ple #1 (where c?pb = dAg «32A). We note that the values 
for \,ujd and /i* are not adjusted, but accepted values 
and so the agreement is quite impressive. However, when 
we used the same parameters for the other samples with 
thinner base layers, we get substantial deviations from 
the model for the first few layers (figure 2). In fact, in 
sample ^1 (and in the theory) the oscillations are super- 
imposed on a global decrease of Tc, in sample #2 the 
background is roughly constant and in sample #3 Tc 
globally increases with the number of bilayers. We argue 
that this behavior is an extension of the inverse proxim- 
ity effect observed in the ultrathin bilayers |0| . The large 
Coulomb interactions characteristics of these thin sam- 
ples are continuously weakened with increasing numbers 
of Pb/Ag bilayers, thus, as more layers are added, we 
approach the Cooper limit proximity effect theory. 

In summary, we have shown that the dc Gcnncs consid- 
erations regarding the Cooper limit are valid for the case 
of multilayers until the layers become very thin. Using 
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FIG. 4: dl/dV-V, normalized to the normal state conduc- 
tance, for sample #1 taken at T=1.65K for each Pb and Ag 
sub-layer. The oscillations of the gap with thickness arc easily 
observed. The inset shows as a function of the number 

kTc 

of layers. Note that this ratio is nearly constant for the last 
few layers. 



[10] The relatively broad transition is a consequence of super- 
conductor fluctuations which introduce some uncertainty 
when evaluating Tc- This can be overcome by fitting 
the curves using the Aslamazov Larkin fluctuation the- 
ory (L.G. Aslamasov and A.I. Larkin, Phys. Lett. 26A, 
238 (1968)). Such fits yield Tcs which are very close to 
the midpoint of the transition as shown in ref 6. 

[11] J.M. Graybeal and M.R. Beasley, Phys. Rev. B29, 4167 
(1984). 

[12] A.E. White, R.C. Dynes and J.P. Garno, Phys. Rev. 

B33, 3549 (1986). 
[13] A.M. Goldman and N. Markovic, Physics Today 51, 39 

(1998) and references within. 
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B40, 6680 (1989). 
[15] P.A. Lee and T.V. Ramakrishnan, Rev. Mod. Phys. 57, 

287 (1985). 



a simple model we were able to deseribe very aecurately 
the variation of the critical temperature and the energy 
gap as a function of the number of layers. We showed 
that for very thin layers the model breaks dow^n indi- 
cating that other physics is involved. We suggest that 
this breakdown is due to the effect of strong Coulomb in- 
teractions in the superconducting ultrathin films. These 
interactions are screened out as normal metal is added to 
the superconductor resulting in an inverse proximity ef- 
fect which must also be taken into account in the thinnest 
layers. 
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